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ABSTRACT 


Nanometer scale fabrication and experimental investigations into the magnetic 
properties of mesoscopic molecular clusters have specifically addressed the need for 
theoretical models to ascertain thermodynamic properties. Technological applications 
germane to these inquiries potentially include minimum scale ferromagnetic data storage 
and quantum computing. The one-dimensional nearest neighbor Heisenberg spin system 
accurately models the energy exchange of certain planar rings of magnetic ions. Seeking 
the partition function from which a host of thermodynamic quantities may be obtained, this 
thesis contrasts two transfer matrix formulations of a classical Heisenberg ring in a 
magnetic field. Following a discussion of the transfer matrix technique in an Ising model 
and a review of material magnetic characteristics, a Heisenberg Hamiltonian development 
establishes the salient integral eigenvalue equation. The 1975 technique of Blume ef al 
turns the integral equation into a matrix eigenvalue equation using Gaussian numerical 
integration. This thesis alternatively proposes an exactly formulated matrix eigenvalue 
equation, deriving the matrix elements by expanding the eigenvectors in a basis of the 
spherical harmonics. Representing the energy coupling of the ring to a magnetic field with 
symmetric or asymmetric transfer operators develops pragmatically distinctive matrix 
elements; the asymmetric yielding a simpler expression. Complete evaluation will require 


follow-on numerical analysis. 
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II. THE PURPOSE OF MAGNETIC MODELS 


A. INTRODUCTION 


In matter, magnetic phenomena both onginate from fundamental forces, and serve 
in ever increasing technological applications. Surprisingly however, there is no single 
comprehensive theory that can fully explain or accurately predict the full variety of 


material magnetic manifestations. For example, the entire theoretical edifice of Quantum 


and Statistical Mechanics cannot answer with certainty detailed questions of 
ferromagnetic coupling. It is known that ferro- and antiferromagnetism arise from short 
range interaction energy that forces spins of unpaired electrons into spontaneous 
alignment. [Ref. 1] More fundamentally, permanent magnetic moments are located in 
atoms or molecules and originate from the circling of the electrons around the nuclei 
(orbital moments) and from the spin of the electrons themselves (spin moments). These 
magnetic moments are proportional to angular momentum which is quantized as integral 


or half integral multiples of 4, (Planck’s constant h~+ 27 ). Particularly in the solid state, 


where electron orbits may be “quenched”, the often intractable many body problem of 
“quantized gyroscope coupling”, or energy exchanges in crystalline lattice structures has 


motivated a variety of models to explain experiment.[Ref. 2] 
1. Historical Snapshot 


In 1907 following Curie’s work, Pierre Weiss [Ref. 3] proposed a theory of 
ferromagnetism in which magnetic moments interact with each other through an artificial 
molecular field proportional to the average magnetization. This type of theory is referred 


to as a “mean-field” theory; mean field theories have only limited accuracy but are often 


useful as a first approach. Subsequent theories incorporated pairwise interaction of 
magnetic moments localized on fixed lattice sites with an energy that achieves a 
maximum value, J, when the moments are either aligned or anti-aligned. Two particular 
models characterizing this energy interaction are the Ising and the Heisenberg models. 
The Ising model assumes the magnetic moments are classical, one dimensional “sticks” 
capable of only two orientations. This mimics the behavior of S=3 quantum spins. Later it 
was found that the Ising model could be applied to a wide range of systems that have an 
essential two-valued nature, such as binary alloys. The Heisenberg model regards the 
magnetic moments as being related to three-component quantum mechanical spin 
operators and assumes the interaction energy is proportional to the scalar product of these 
operators. There are other “spin” dimensional models but only a select few have been 
solved exactly for various space dimensionalities and external magnetic fields. Stanley 
{Ref. 4] provides an excellent comprehensive summary of these models, their 
applicability and limitations. 

Neither the Ising nor the Heisenberg model has yielded as yet to an exact solution 
for a three dimensional (3-space) lattice. In 1944, in a landmark in the history of phase 
transitions and critical phenomena, Onsager [Ref. 5] solved the two-dimensional Ising 
model. An infinite-spin version of the 1-D Heisenberg model with free boundary 
conditions was solved by Fisher [Ref. 6] in 1964. Fisher showed that an infinite-spin 


Heisenberg model was equivalent to a classical version of the Heisenberg model in which 
the quantum spin operators are replaced by classical vectors of length ./S (S +1) that are 


free to orient in any direction. This classical counterpart to the quantum Heisenberg 


model is called the classical Heisenberg model. The classical Heisenberg model should 


apply to high-spin magnetic systems for all but extremely low temperatures. Fisher’s 
student Joyce [Ref. 7] in 1967 published an exact solution to the zero-field one 
dimensional isotropic classical Heisenberg model with cyclic boundary conditions 
employing Wigner 37 symbols. Blume et al [Ref. 8] in 1975 employed a transfer-matrix 
integral equation method and extended Joyce’s work to tackle a 1-D classical Heisenberg 
system in an applied magnetic field. This thesis will follow and compare this latter 


method with a currently proposed method by Auslender [Ref. 9]. 


2. Models and Thermodynamics 


These Ising and Heisenberg “toy models” enable reasonably accurate theoretical 
descriptions of certain physical systems and, significantly, shed valuable thermodynamic 
insight on some fluid and magnetic phase transformations. Statistical mechanics 
establishes a connection between the microscopic and macroscopic, or bulk, 
thermodynamic descriptions of a system. Central to statistical mechanical formulations is 


the partition function, 





Jb = > exp(- 24¢(c)) , where B= (1.1) 


kg_ is the Boltzmann constant, 7 is the absolute temperature, and HH is the energy 


Hamiltonian for each available quantum state, o. As a summation of all Boltzmann 
factors, the partition function 1s the inverse proportionality factor between the probability 


of a particular energy state, P(o) and each Boltzmann factor, 


eo felc) 
Z 





P(c)= (le2) 


If ascertainable, the partition function is a very useful result. In fact, the partition 
function can be called the holy grail of equilibrium statistical mechanics because 
essentially an entire thermodynamic description of a system can be derived from this 
function. The generalized ensemble theory of Gibbs enables computing the complete set 
of thermodynamic quantities fom purely mechanical properties of its microscopic 
constituents assuming only a “mechanical” structure, and obedience of Lagrange and 
Hamilton’s equations of motion [Ref. 10]. Thermodynamic averages such as entropy, 
average energy, heat capacity, magnetization and susceptibility, as well as the Gibbs 
potential and particularly the Helmholtz free energy are derived directly from the 
partition function. Of course, the essential completeness of the partition function 
necessitates a summation over all states; therefore, obtaining the partition function is no 
small challenge. This underscores the value of models that permit an exact determination 
of the partition function. A model can enable a theoretical calculation of state variables. 


Of course, a model’s validity and utility corresponds directly with its fidelity to reality. 


B. MESOSCOPIC MOLECULAR CLUSTERS 


Recently, fascinating experiments with large molecular clusters of metal ions 
provide an opportunity to employ a one dimensional Heisenberg spin model. Of interest 
is both the extremely clever and revealing experimental techniques and the particular 
scale (nanometer) of these investigations, where renormalization group theory [Ref. 11] 
had heretofore “coarse grained”, between atomic and bulk scales. These “mesoscopic” 
magnetic molecular clusters are enabling investigation of such behavior as quantum 
tunnelling of magnetization [Refs. 12, 13]. Some authors have forecast technological 


applications of these ultrasmall complexes in the field of both data storage and quantum 


computing [Refs. 14, 15, 16] . Two molecular structures, depicted in Figure 1, are 
particularly noteworthy for their symmetry, high spin and revealing characteristics. The 
first contains twelve manganese ions, arranged in a ring of eight Mn” ions with spin S=2 
aligned parallel, enveloping the remaining four manganese Mn* atoms which form a 


tetrahedron with spins S=3 in the opposite direction to the encircling octagon. This 


manganese acetate, Mnj2 O12(CH3 COO)16 (H2 O)4 is described as superparamagnetic, 





Schematic view of the core of a 
[Mn,.0,.(carboxyiato), .} cluster in which only the 
metal atoms and tne bridging oxygen atoms 
(small circles) are shown. The manganese(IV) at- J 
ams are enhanced by the shadowing. _ - 


View of the ring structure of the Fe,, Clus- 
ter, where the dotted circles represent the iron 
atoms and the empty circles are, in order of de- 
creasing size, chlorine, oxygen, ang casoon. 


Figure 1. Two Mesoscopic molecular clusters, the manganese acetate on the left and the 
ferric wheel, right. (From Ref. 15) 


having a ground state of S=10 and the measured spin dynamics of this cluster in varying 
magnetic field reveals a hysteretic magnetic relaxation ascribed to resonant tunneling 
between quantum spin states. [Refs. 17, 18, 19] The second noteworthy molecular 
cluster is [Fe(OCH3)2(O02CCH2CL)]i0, known as a “ferric wheel”. It contains 10 nearly 


coplanar Fe*” ions each of spin S= 5, symetrically positioned on what constitutes a 


Heisenberg ring, a planar, one-dimensional spin system. In addition to this “ferric wheel” 
decagon, other iron molecular clusters have been synthesized, (e.g. Feg , Fei7 ,Fio, as 


well as Fejo). 


These molecular clusters would therefore seem to contain sufficiently few 
magnetic constituents that an exact determination of the partition function could be 
attempted. According to Gatteschi et al, however, [Ref. 15] when commenting on the 
quantum mechanical energy computations of these ferric rings: 


A quantitative interpretation of the magnetic properties of these 
compounds has been possible only for Feg, and this at the cost of some 
effort [Ref. 20], the total number of states being 1,679,616. Exploiting 
symmetry allows the reduction of the problem to that of calculating 81 
matrices, ranging in dimension from 1 to 4,170. A similar analysis proved 
to be impossible for the Fe)7 or Fejs clusters. Even with use of all the 
possible symmetnes, the dimensions of the matrices remain much too 
large to be tackled with the standard approach. In this field, theoretical 
developments are strongly needed, so that we can interpret the 
thermodynamic properties of the new materials. 


This clarion invitation to employ the theoretical Heisenberg spin model has been 
answered by Luscombe ef al, [Ref. 21] astutely developing approximations for the 
relevant thermodynamic qualities. In [Ref. 21] 1t was shown that the classical Heisenberg 


model well approximated the observed thermodynamic properties of the “ferric wheel” 


cluster with 10 Fe®* ions of spin S=5/2. This collaboration between theoretician and 
experimentalist was acknowledged by the latter in January, 1998, when Lasciafari et al, 
[Ref. 22] the leading research team in the field, effectively employed theoretical 
thermodynamic results from [Ref. 21] to advance the macromolecular magnetic frontier. 
The continued development of these molecular magnetic systems, hopefully into 
useful nanomagnetic technologies, will clearly require reliable and robust theoretical 
techniques for predicting their thermodynamic properties. As noted above, it has been 
shown [Refs. 21 & 22] that approximate treatments based on classical Heisenberg spins 
can predict extremely well the observed magnetic behavior of small quantum Heisenberg 
systems. It is thus worthwhile to search for improved theoretical approaches to modeling 
the thermodynamic properties of classical Heisenberg spin systems. As will be discussed 
below, Blume’s method [Ref. 8] for evaluating the partition function of the 1-D classical 
Heisenberg model in an applied magnetic field involves solving numerically an 
eigenvalue integral equation based on what 1s known as the transfer-matrix operator as its 
kernel. Using a numerical Gaussian integration technique, Blume ef al transform the 
eigenvalue integral equation into a matrix eigenvalue equation. Auslender, [Ref. 9], has 
recently proposed an alternate strategy for solving the integral equation. Auslender’s 
proposal is to represent the transfer-matrix operator in a basis set of spherical harmonics. 
As will be shown below, the spherical harmonics are the eigenvectors of the transfer- 
matrix operator for zero magnetic field in this basis set. Whether Auslender’s proposal 
results in a more efficient method from a numerical point of view remains to be seen. In 


this thesis, we will, for the first time, set up the matrix representation of the transfer- 


matrix operator in the spherical harmonic basis set. Since the transfer matrix method is 


key to the results of this thesis, we will review this method in the next section. 


C TRANSFER MATRIX FOR AN ISING MODEL 

Attributed to Kramers and Wannier [Ref. 23], the simplest illustration of the 
transfer-matrix technique is its application to the one dimensional nearest neighbor N- 
spin Ising model in a magnetic field. The physical justification will be presented in the 
next chapter, what follows simply demonstrates the math. 


The 1-D nearest neighbor Ising Hamiltonian is defined by 


N H N 
KH =—-J PF Fi41 ser +0,4), ( 1.3) 
i=] i=] 


where o, = +1 isthe randomly up or down onented spin 3 Ising variable at lattice site i, 


1<i<N, and where J is the nearest neighbor exchange parameter, and H is the applied 
static magnetic field which is parallel or antiparallel to the moment of each “invertable” 
spin. Although an open Ising chain is solvable without recourse to the transfer-matrix 
method (e.g.Stanley [Ref. 4]), here we will instead assume periodic boundary conditions 


and define. o,,,, =0,. Recalling equation (1.2), thermodynamic averages are constructed 
from the probability distribution P(o)=Z, exp(— A2¢(c)), where 8 =(k,7)., 


and Zy 1s the partition function, 


Zy(K,L)=¥ exo Bron o.))= Dew KL + FDC, +e.) (4 
{co} {co} i=l 


i=] 
The notation {co} indicates a summation over 2" spin configuration, i.e. 


1 


1 
=> 2 >, , and K = J and L = 6H are dimensionless coupling constants. 


{o} 3, =-1 oy=-l 


Using the fact that Ising variables are classical “sticks” that will commute, the 


exponential of the sum ts a product of the exponentials and (1.4) may be written, 


Z (Kee > T(e,0,)T(62,03) a TCy.1,6y T(o,2,), v5) 
{o} 


ib ’ hae 
where 7 (c,0,,,) = exp[Ko,a,,, + ra + o,,,)] which written out explicitly is a 2x2 


transfer matrix with elements T(c,o’), each spin o having possible values + 1, 


T(+14+1) T(+1-1)) (e%* e* 
fe = = - : (1.6) 
T(-1+1) TG1-1)) le* ek 
Then, summing over spins 2,...,N in (1.5) the partition function is given by 
} 
ipigpe > (7% \o,0,) = Trace en (1.7) 


o,;=-1 
that is the partition function of the N-spin Ising model with periodic boundary conditions 


is given as the trace of the N™ power of the transfer matrix. Since the trace of a matrix is 
equal to the sum of its eigenvalues, and the eigenvalues of T‘ are A” and 2% where these 


are determined by the equation, 


e e- ~-—jA 





} 
= 0 with solutions a =e* cosh L+(e* sinh? L +e )2 (1.8) 





we finally arrive at the result Zy= 7," +2" = 2 1 (2 (1.9) 


Since 41>A2, the second term in the parenthesis in (1.9) goes to zero for large N and can 


be neglected. In zero magnetic field, H=0=L, (1.9) including both eigenvalues, yields 
Zy = 2" (Cosh"K + Sinh’K). (1.10) 
This nearest neighbor Ising discrete spin system engenders a “2 by 2” transfer 
matrix. The nearest neighbor classical Heisenberg model however, has a continuously 
directable, all aspect spin system and the transfer “matrix” for continuous spins is an 
infinite matrix — or the kernel of an integral equation. Nevertheless, the eigenvalues of 
the transfer matrix allow one to obtain the partition function. We will discuss in Chapter 


IJ the application of the transfer-matrix method to the classical Heisenberg model. 


D. THESIS OBJECTIVES 

This thesis will contrast two methods of solving the one dimensional classical 
Heisenberg spin model. Blume and Auslender are the authors of the two techniques, the 
former acknowledged as the accepted method and the latter, a proposed alternative. Both 
Auslender’s and Blume’s methods are concerned with obtaining eigenvalues to the 
transfer matrix associated with the one-dimensional, nearest neighbor classical 
Heisenberg model in a magnetic field. Both methods seek the solutions of the eigenvalue 


equation: 


fas’ T(s,5') Vi n(S')= Amin (S); (1.11) 
where 7 (s, s’) = ex Ks -s'+ = Us, ae )} ds = sin @d@d¢, and s,s’ are adjacent three 


component spin vectors. Both methods seek the partition function for the N- spin system 


with periodic boundary conditions, which is given in terms of the eigenvalues 


Zy => Aim * (1.12) 
l,m 
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For zero magnetic field, the eigenvalues and eigenvectors are exactly obtainable as first 


shown by Joyce. [Ref. 7] In this case 4, ,, = 427, (K ) where /; (K) are the modified 
spherical Bessel functions, and y,,, Behe (0,9), the spherical harmonics. 


Blume et al, (Heller and Lurie coauthors)[Ref.8], turn the integral eigenvalue 
equation into an M x M matrix eigenvalue equation using M-point Gaussian integration. 
With a matrix size of 16x16, this technique results in a convergence to seven significant 
figures for the values of the Hamiltonian. Auslender also turns the integral equation into 
a matrix equation. The matrix in this case is obtained by expanding the eigenvectors 
using the spherical harmonics as a basis set. That is, he suggests expanding the 
eigenvectors in a magnetic field using zero-field eigenvectors. In principle, Auslender’s 
matrix is infinite dimensional and must be truncated at some point. It is likely that a 
comparison of merit will only result from detailed numerical computations; however, this 
thesis will merely discuss formulations that reduce the integral equations to matrix 
eigenvalue equations employing both methods. Leading up to these formulations, any 
discussion of equilibrium statistical mechanics of the classical Heisenberg spin model 
must commence with formulating the energy exchange Hamiltonian. The transfer matrix 
will then lead to setting up the eigenvalue equation. There will be some discussion of 
zero magnetic field behavior and necessary discussion of integral equations, especially 
the role of a symmetric kernel. Representing the energy coupling of the Heisenberg ring 
to a magnetic field with symmetric or asymmetric transfer operators will be shown to 


result in mathematically equivalent but pragmatically distinctive element formulations. 


1] 


Prior to this theoretical development in Chapter IIL Chapter I will consist ofa 
cursory review of the magnetic properties of matter following largely the excellent if old 
treatments of Von Hippel [Refs. 1&24], Kittel [Ref. 25], and Ashcroft & Mermin [Ref. 


26]. 


| 7a 


U. MAGNETIC PROPERTIES OF SOLIDS 

A. OVERVIEW 

The magnetic properties of solids onginate in the motion of the electrons and in 
the permanent magnetic moments of the atoms and electrons. This chapter will provide a 
review of magnetic characteristics and hopefully lead to a “motivation” of the Heisenberg 
model. Diamagnetism, which is very weak, arises from changes 1n the atomic orbital 
states induced by an applied magnetic field. Paramagnetism results from the presence of 
permanent atomic or electronic magnetic moments. Ferromagnetism, which is very 
strong, occurs when quantum mechanical exchange interactions align adjacent magnetic 
moments in the same direction. Ifthe erehange interaction aligns the moments in 
opposite directions, and only one type of moment is present, cancellation occurs and the 
material is called anti-ferromagnetic. If two or more types of moments are present, there 
is a net moment equal to the difference and the material 1s called ferrimagnetic. Above 
some Critical temperature, a phase transformation occurs and a ferro-, antiferro-, or 
ferrimagnetic material becomes paramagnetic. Ferro- and ferrimagnetic materials consist 
of domains or regions of completely magnetized material, separated by boundaries 
known as Bloch walls. According to Kittel [Ref. 25], and Hippel [Ref. 1], domain 
structure, dynamics, and boundary displacements are determined by various types of 
energies, such as magnetostatic energy, crystal anisotropy and magnetorestrictive energy. 
The complexity of these resultant forces contribute to the scientific and technological 
richness in this field. The succeeding paragraphs merely scratch the surface of these 


topics. 


B. ANGULAR MOMENTUM AND MAGNETISM 
The relation between angular momentum and magnetism is based on the mac- 
roscopic observation that a current / circling an area A creates a magnetic field identical 


to that of a magnetic dipole. As such, the Bohr hydrogen atom’s magnetic dipole, 


|u| = JA = evar’, (2.1) 
applies for an electron circling the proton v - times per second in an orbit of radius 7. 
For this same orbit the classical mechanical angular momentum, 

IL] = |v xr| = m2arv r. (2.2) 


This angular momentum L is antiparallel to 44 and combining (2.1) and (2.2) shows the 


; e 
magnetic and angular moments are related as w= ar , 
m 


Thus the magnetic and mechanical moments of circling electrons are interdependent and 


the gyromagnetic ratio y is classically defined y = fe ~ om , (2.3) 
m 


At atomic scales, the Bohr magneton 1s considered an elementary magnetic moment with 


| Ls ele vi =9.27X10-* in units amperes meter’ or joules/tesla (2.4) 


2m, 


(Note: of course, in measuring magnetic moments of nuclei, the nuclear magneton would 
: eee l 
be a preferred unit and with a mass substitution, | Ly| = 7336 He l) 


If the magnetic moment is measured in Bohr magnetons and angular momentum in units 


of f, the ratio of magnetic to mechanical moment, known as the dimensionless g-factor, 


a h h 
(for the classically single orbiting electron), g =——— =y—_ = 1 
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In a magnetic field B, the permanent magnetic moments will experience a torque 


een o, but p= eee eS BS yee (2.5) 
at e dt 2m 


For a static field applied in the +z direction, taking the cross product followed by the 
second time derivatives, yields the component equations of (2.5): 

ii, =-(B) un, . ii, =-(BYu,,  #, =0. (2.6) 
Solutions of (2.6) are: 4, = Acosw,t, bw, =ASinw,t, and pu, =const, where 

@, =yB andw,/2a =v, is called the Larmor frequency. Hence these two oscillating 


components are 90° out of phase and add to a circular rotation in the x-y plane. The 
magnetic moment precesses around the magnetic field axis with a frequency proportional 
to the field strength but independent of position. 


The Larmor frequency is not quantized but three quantum numbers are integer 


multiples of i. The boundary conditions on the time independent Schrodinger wave 


equation restrict the quantum numbers as follows. The principle quantum number 7 is 


allowed positive integer values 1,2,3,... The orbital angular momentum quantum number 


é can take integer values 0S 7<n. Quantum mechanically the total orbital angular 


momentum =L=, e(@+1\i. The magnetic field directed component of L=l,=mh, 


where the magnetic quantum number, 


m=./£(@+1) cosé@, (Zi) 


with M= £,2-1).0..0.... —(£ ~ 1), —£,and @ the quantized polar angle between the B field 


and L. 
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The quantized magnetic moment in the magnetic field has a potential energy 
=—y-B= pBcosd =—Je(e +1 ity |B|Cos@ =my,B. (2.8) 
The electron itself, has intrinsic angular momentum and thus creates a magnetic 


rs 4 39 


e . 
moment W =~-g, =e where the electronic“g.” factor was predicted by Dirac to equal 
m 


approximately 2, has been measured experimentally to 2.0023, and is given [Ref.26] by, 


2 
a (2.9) 
4né,hc 137 





Z,= a1 + = + O(a? )+ Bare | , a, the fine structure constant = 
a 


| ee Pd 


Using the electronic spin “g” factor equal to twice the orbital, then classically the 


permanent magnetic dipole moment yw =—y(L+2S). (2.10) 


Cc. SPIN ORBIT COUPLING, SHELL MODEL, AND HUND’S RULES 


The combined angular momentum produced by the spin and orbital motion is 
J=L+S. The total angular momentum J 1s always a good quantum number, (i.e. 
commutes with the Hamiltonian), but L and S are good only to the extent that spin — orbit 
coupling is unimportant[Ref. 26]. Both the spin and orbital angular moments will tend to 


precess around J and, 


J = (L+S) = 245° +2L-S=e(e +1 +20 +2L-S (2.11) 


The vector addition first solved by Lande results in the g — factor for Z-S coupling 


J(J +1) + S(S +1)- L(L +1) 


=(L+2S)eJ/J* =1 
eter ° 2 (J +1) 


(2.12) 


where this Landé g- factor is relevant in the expression u=-gu, J/h. 
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Glibly allowing that the colossal variety of not only magnetic effects, but of all 
nature’s splendor derives quantum mechanically from the atomic shell model, angular 
momenta coupling and shell filling configurations are key to understanding this process. 
Deriving from the Pauli exclusion principle or the antisymmetry of fermion wave 
functions, the underlying quantum mechanical justification of shell filling, selection 
rules, ionization potentials, electron affinities and atomic bonds are beyond the scope of 
this thesis, but brief essentials relevant to magnetic properties follow. Filled shells will 
have zero orbital, spin, and total angular momentum, and consequently zero permanent 
magnetic moment. For ground states in partially filled shells, the order of filling obey the 
Pauli exclusion principle and is roughly agains: by Hund’s rules which are as follows: 

Rule 1. In silasine n electrons into the 2(2£ +1) levels of the partially filled shell, 

those that lie lowest in energy have the largest total spin S; thus if possible, the 

first 22+1 of allowed electrons in a shell will align spin-up. 

Rule 2. The total orbital angular momentum L of the lowest lying states has the 


largest value that is consistent with Rule 1 and the exclusion principle 


Rule 3. Total angular momentum J = JJ (J +1) where J takes on integral values 
between |Z — S| if the shell is less than half full and L +S if more than half full. 


Russel-Saunders coupling applicable to lighter elements and favored in the d-shell and f- 
shell, is represented by a term in the Hamiltonian of the form A(L-S) This spin-orbit 
coupling will favor maximum J (parallel orbital and spin angular momenta) if A is 
negative, and minimum J (antiparallel orbital and spin momenta) if 4 is positive. As it 


turns out, 2 is positive for shell that are less than half filled and negative for shells more 


iy 


than half filled. The configurations of ground state d-shell and f-shell ions are tabulated 
below. The “multiplets” notation "***1X,", where (X=L in the ‘SPDF’ spectroscopic 


code), conveys S, L, and J, and the shell filling order for solids of magnetic interest 


illustrate Hund’s rules. 


d-shell (I = 2) 
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Table 1. Ground states of ions with partial d- or f-shells per Hund’s rules* [From Ref.26] 
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Bb. MAGNETIZATION AND DIAMAGNETISM 


The magnetic induction B in free space is related to the field strength or magnetic 
intensity by B= ,H where yu, =42x107’ henry/meter is called the permeability of 
vacuum. Ina solid material B = wH, which alternatively can be expressed as, 

B= y,(H+M)= ue. (2.13) 
(Note: Of course, u here is the permeability of the solid not a dipole moment.) M is 
called the magnetization of the solid, the term 4,.M equals the extra magnetic induction 
due to the material. M in fact, is equivalent to the density of magnetic dipole moment or 
dipole moment per unit volume. The magnetization is also proportional to the applied 
field and the factor of proportionality is called the susceptibility. The magnetic 
susceptibility per unit volume is defined as y=M/H. (2.14) 

Substances with a negative magnetic susceptibility demonstrate diamagnetism 
which is a material manifestation of Lenz’s Law, which in effect orients Faraday 
induction such that “a current induced by a changing field will always oppose the change 
that induces it.” With zero angular momentum, fully closed shells have zero permanent 
moments, (eq 2.10), but in an external magnetic field, there is an induced moment. The 


averaged induced magnetic moment for each electron equals, 


e?(r\’ 
HB induced = ae ? (2.15) 
m 


where e and m equal the electron charge and mass and (7,)1s the average electron orbital 
radius. Diamagnetism in most solids is very weak with susceptibilities on the order of 
10°. Itis generally only observed when other types of magnetism are totally absent. The 


variety of net magnetic susceptibilities in matter is illustrated by Figure 2. 
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Figure 2. Diamagnetic elements in lower portion, paramagnetic in upper [From Ref. 27] 


: PARAMAGNETISM 


Positive susceptibilities (y>0) are generally termed paramagnetic. Permanent 


magnetic moments give rise to paramagnetism, and per Kittel [Ref. 25], electronic 
paramagnetism is found in: 
a) All atoms and molecules possessing an odd number of electrons, since the 
total spin of the system cannot go to zero. 
b) All free atoms and ions with a partly filled inner shell: transition elements, 
rare earth and actinide elements. 
c) A few miscellaneous compounds with an even number of electrons, including 
molecular oxygen and organic biradicals. 


d) Most but not all metals as depicted in Figure 2. 


What follows is a Maxwell Boltzmann distribution treatment of the Langevin Theory of 
paramagnetism. 

1. Langevin Function 

Permanent magnetic moments tend to orient in magnetic fields. With N atoms per 
unit volume, each bearing a magnetic moment 44 magnetization results from the 
Orientation of these moments in an applied field. Thermal disorder resists this orientation 
tendency. The energy of interaction with an applied magnetic field His 


E=-u-H =—pAHcos@, where @ is between the moment and the field direction. The 


magnetization will be M=Ny cos@ where N is the density and cos@ is the average over 
a distribution in thermal equilibrium. According to Boltzmann distribution, the relative 


probability of finding a molecule in a solid angle element @2 is proportional to 


oe , and cos@ = fe *cos@iQ +(e *aQ Over all solid angles, 
cos0 = | 2zsin6 cos eM! 9M 1g =. [2x sin A etties Ol 1g (2.16) 
0 0 


letting x=cos@ and a= yA /kT, then 


] 1 d 1 l 
9 = |e*xdx + |e" de = —In | e“dx = cotha—-—= Lia). 2.17 
cos fe fe a n fe cotha - (a) (2.17) 


-1 -1 “1 
L(a) is called the Langevin function. When the field energy is small in comparison with 
kT, a<<1, then L(a)~a/3 and M~»x Nu?H/3kT (2.18) 
The magnetic susceptibility in the limit as wH/kT<<J 1s 

vy =M/H =Nu’? /3kT =C/T, (2.19) 


where the Curie constant C = Nu? /3k. The inverse temperature dependence is known 
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as the Curie law and the entire expression is called the Langevin equation. This Langevin 
derivation is entirely classical with unrestricted space orientation of the moments in a 
magnetic field and furthermore, depends intrinsically on the Maxwell Boltzmann 
distribution. A quantum theory of paramagnetism still employing the Boltzmann 


distribution uses the Landé g factor (2.12) and what is known as the Brillouin function 
for calculating the 2/+/ discrete and equally placed energy levels in the field. 


Essentially equivalent to the Curie Law, the calculation yields y = Np* yu," / 3kT , where 


the effective number of Bohr magnetons is defined as p = gl (U + 1p 

2. Pauli Paramagnetism 

The Langevin equation does not apply to conduction electrons which obey the 
Fermi-Dirac distribution. Conduction electrons are neither spatially localized like 
electrons in partially filled ionic shells, nor because of stringent constraints of the 
exclusion principle, do they respond independently like electrons localized on different 
ions [Ref. 26]. Although small, Pauli paramagnetic susceptibility results from the 
coupling of intrinsic electron spins with an applied field H. There is also a diamagnetic 
effect arising from the coupling of the field to the orbital electron motion. This is called 


Landau diamagnetism and for free electrons in metals, the susceptibility, 


l pee : 
V7 onda a XPauii he resulting net susceptibility for V conduction electron is 


y=Nu, / FE’, , where Er is the Fermi level. Pauli paramagnetism is independent of 


temperature and even at room temperature 1s hundreds of times smaller than the 
paramagnetism of magnetic ions. Paramagnetism usually masks the atomic 


diamagnetism present in solids. In practice, it is the total susceptibility that is revealed by 
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a measurement of bulk moment induced by a field and this is a combination of the Pauli 
paramagnetic susceptibility, the Landau diamagnetic susceptibility, and the Larmor 
diamagnetic susceptibility (of the closed-shell ion cores). As a result, isolating 
experimentally these particular terms of the susceptibility is not at all straightforward. 
Nuclear magnetic resonance (NMR) is one such technique that enables experimental 
discrimination of these susceptibilities. Like NMR which can measure spin-lattice 
relaxation rates, another technique called muon spin relaxation, also is central to the 
current investigations of magnetic molecular clusters cited in the introduction [Ref. 22]. 
Unlike these recent frontier if somewhat esoteric inquiries, the next section will attempt 


to describe a more prosaic phenomena, namely refrigerator magnets. 


| FERROMAGNETISM 


The transition metals Fe, Co and Ni, rare earth metals such as Gd and a few 
oxides such as CrO2 and ErO display very large magnetization. These ferromagnetic 
materials contain permanent atomic magnetic dipoles, the difference from a paramagnetic 
substance being that, below a certain temperature, the dipoles retain parallel orientation 
even in the absence of an external field. Figure 3a. depicts a magnetization curve ofa 
ferromagnetic material. This hysteresis loop characterizes the magnetic induction B as a 
“function” of the applied field H. As the applied field H is increased, B begins to 
increase slowly. The slope rises sharply as B rapidly increases until the saturation 
induction. Upon decreasing the field, the original curve is not retraced. At H equal to 
zero, the specimen is still magnetized with the remanent induction. Here is the reason 


that zero field permanent magnets are able to emblazon refrigerators. If H is now made 
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negative, when B=0 indicates the coercive force required to de-magnetize the material. 
The symmetric curve depicts saturation, remanence and coercive force for negative 
induction values. This irreversible double valued hysteresis is the signature behavior of 
ferromagnetic materials. The work required to go around the hysteresis loop 

once 1s proportional to the enclosed area. Technologically an ii with a fat loop 
(Figure 3b.) makes a good permanent magnet; whereas a thin loop, (Figure 3c.) with 


small area, demagnetizes rapidly and makes an efficient AC transformer element. 
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Fi gure ). Hysteresis curves for (a) soft iron, (b) a good permanent magnet, and (c) an 
alloy suitable for use in a power transformer. [From Ref. 27] 


The source of ferromagnetism is a parallel alignment of unpaired electron spins. 
As noted in the introduction, Weiss (1907) postulated a molecular field to explain ferro- 
magnetism and he further postulated domain formation to explain the hysteretic magnet- 
ization curve. The molecular Weiss field was formulated as Hw=AM, where A is termed 


the Weiss constant and this field is added to the applied field H in Curie’s Law (eq 2.18). 
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Then the magnetization, 


C | IM| mM G C 
ae AM d sol for —- =—=y= = 220 
M = (H+ ), an solving for Hl = Oe TAC T-T.” (2720) 








where Tc the Cune temperature, is where the transformation from the paramagnetic to the 
ferromagnetic occurs. The domain hypothesis can be inferred from the hysteresis curve 
and were observed in fact some twenty five years later in 1931. Within domains, all 
moments are aligned but there is random orientation of the domains resulting in a net 
magnetization of zero. The external magnetic field induces magnetization via domain 
wall motion. The applied external field will cause favorably aligned domains to grow 
thereby shrinking unfavorably oriented domains. Rotation of other domain moments 
maximizes the magnetization. The saturation value corresponds essentially to single 
domain status. When the field 1s removed, the specimen remains magnetized. Although 
domains typically tend to rotate back, the large aligned domains do not easily revert to 
the original random arrangement. Reduction and reversal of the field allow a domain 
pattern to return, depending on the ease with which domain walls can nucleate, move 


through the material, and again be ejected. 


G. ANTIFERROMAGNETISM AND FERRIMAGNETISM 


There are two other important classes of magnetic behavior. When adjacent 
unpaired spins are aligned in an opposite fashion, the resultant phenomenon is called 
antiferromagnetism. The susceptibility 1s then positive and increases as the temperature 
increases since thermal energy as always disrupts long range order. Figure 4. compares 
alignment of magnetic moments and temperature effects on ferro- and antiferromagnetic 


materials. The peak in the susceptibility y versus temperature T in Figure 4b. is called 


os 


the Neel temperature Oy and corresponds to the Cune temperature in ferromagnetic 
materials. A class of Rit oat oxides known as spinels, having the composition 
XOFe203 (where X is a metal), exhibit ferromagnetic interaction yet have anti-parallel 
spins as depicted in Figure 4c. A net moment results since the opposite spins are 


unequal. The magnetization of these spinels, known as ferrites, have wide application in 


the electronics industry. 
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Figure 4. Magnetic susceptibility versus T for (a) ferromagnetic, (b) antiferro- 
Magnetic, and (c) ferrimagnetic materials, with magnetic moment alignments 
indicated for each case. [After Ref. 1] 


H. QUENCHED ORBITS AND THE HEISENBERG MODEL 


The moments of ferromagnetic arrays could in principle stem from orbital 
moments as well as spin moments of individual electrons. However, since crystal 
structures are held together by electron bonds, it is not surprising to discover from 
magnetomechanical measurements that the orbital moments are essentially quenched by 


such bond formation. The gyromagnetic ratio for orbital moments 1s e/2m (eq. 2.3); for 
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ferromagnetics, it proves to be about e/m, the value of electron spins. Although the 
contribution of orbital moments to the saturation magnetization is generally only 5-10% 
of the total, they are roneneless important as a source of magnetic crystal anisotropy 
[Ref. 1]. A quantum mechanical explanation for the Weiss field proposed by Heisenberg 
(1928), involves an exchange interaction between neighboring electron spins. The 
exchange energy explanation may be motivated by the Pauli exclusion principle and 
corresponding Fermi Dirac statistics requiring distinction af each electronic state by its 
own unique set of quantum numbers. Overlapping wave functions can lead to a decrease 
in over-all energy in certain cases, and therefore favor a parallel alignment of spins. The 
spin quantum number corresponds to up or down, hence inversion from parallel to 
antiparallel leads to a new electron cloud of different electrostatic energy. Bohm 
[Ref.28] concurs and attributes the antisymmetry of the complete. electronic wave 
function with prescribing parallel or antiparallel spin alignments. Furthermore the 
energy, apparently a result of spin interactions, 1s actually a result of the correlation 
between mean coulomb energy and spin. 


Ashcroft and Mermin [Ref. 26] construct a spin Hamiltonian for a two electron 
, ae | , 1(1 3 
system noting that each individual electron spin operator satisfies S, = B ae = Y, , 
so that total S satisfies , 
3 
S? = (Seas =5+2S, -S;, (221) 
since S* has eigenvalues S(S+1) in states of spin S, the operator S, -S, has eigenvalue # 


in the singlet (S=O) state and +3 in the triplet (S=1) states. Consequently the operator 


Zi 


wera“ (E, +3E,)-(E, -E,)S,-S, has eigenvalue Es in the singlet and E, in each of 


the triplet states and is the desired Hamiltonian. By redefining the zero of energy, the 


constant, common to all four states, (E, +3E,)/4 can be omitted and the spin 


Hamiltonian is #7" =-JS,-S,, J=E,-E,. (2.22) 


Seeking lowest energy, the scalar product of the vector spin operators will favor parallel 
spins if J is positive and antiparallel if J is negative. It 1s also noteworthy that in contrast 
to dipolar interaction, the coupling in this spin Hamiltonian depends only on the relative 
onentation of the two spins and not on the vector difference between the spins. It is 
remarkably true that in many cases of interest, the form of the spin Hamiltonian is simply 


that for the two spin case summed over all pairs of ions, 
i, SS, | (2.23) 


This expression (2.23) is called the Heisenberg Hamiltonian and the J, are the exchange 


coupling constants. Stanley [Ref. 4] points out that this model 1s not valid for a wide 


variety of real magnetic materials as it assumes: 


a) Well localized spins (i.e. small wave function overlap) 


b) Complete isotropy of interaction. 


The 3d transition metals have overlapping wave function and rare earth metals are 
generally anisotropic. Nonetheless, the Heisenberg spin Hamiltonian formulation can 
yield fruitful theoretical information in many cases, one such being the mesoscopic 


“Ferric” Wheels. 


28 


Ii. HEISENBERG SPIN SYSTEM FORMULATIONS 


A. OUTLINE 


This chapter will contrast two theoretical approaches to obtaining the partition 
function for a one-dimensional nearest neighbor classical Heisenberg spin system in a 
magnetic field. At the outset, section B will develop the symmetric and non-symmetric 
transfer operators, and show that the partition function for an N-spin system is equal to 
the sum of the “N‘"-power raised” transfer matrix eigenvalues. Next section C will 
delineate the zero field analytic solution of a classical Heisenberg ring. Then, tackling a 
finite magnetic field in Section D, the numerical integration approach of Blume ef al will 
be described, following directly the authors’ formulation. Section E will suggest a new 
approach to solving the classical Heisenberg spin eysten The matrix eigenvalue 
equation will be constructed by representing the transfer operator kernel in a basis of 
spherical harmonics, (which are the zero field eigenvectors in the analytic solution). The 
infinite matrix that results will be examined qualitatively for both zero and non-zero field 
characteristics. Development using first a symmetric magnetic field transfer operator, 
followed by and compared with a non-symmetric transfer operator expansion, employ 
both integral equation and rotation group mathematics. The non-symmetric transfer 
operator surprisingly yields a simpler matrix construction and both fommnierions enable a 
conceptual contrasting to Blume et a/ while setting the stage for an actual calculational 


comparison. 
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B. TRANSFER MATRIX FOR CLASSICAL HEISENBERG SPINS. 


The starting point in this development is the Heisenberg Hamiltonian, which is a 


summation of quantum spin operators in units of i: 


N 
H=—-J » Ss; -S,,,, where S,,, =S, (periodic boundary conditions), and J is the 


i=l 
unique exchange interaction energy applicable to each adjacent pair. It is noteworthy to 
recall equation (2.22), that in this formulation (/<0) J>0 promotes (anti-) ferromagnetic 
ordering at low temperatures. To incorporate system coupling with the magnetic field, 


the potential energy of a magnetic moment yw ina magnetic field B, is -u-B. Quantum 
mechanically, u=-guzJ /h where g is the Landé g-factor, wp is the Bohr magneton and J 


is the total angular momentum. This treatment will assume that orbital angular 


momentum L is completely quenched by the crystal fields (1.e. L=0), so the potential 


; | 
energy term, (also in units of A), is given by guz,B- > S,. Ifwe take the B field as 


i=l 


N 
defining the z-axis, then the magnetic field energy term 1s given by gu ee S? , where 


i=l 
S” is the z-component of the spin at site i. For later convenience and without any loss of 


generality, we define m=-—gy, and thus the total Hamiltonian 1s given by 


N N 
H =-J>S,-8;,,-mB> S; . Gal) 


t=] t=] 
Now the classical spin approximation, (which 1s necessarily invalid at low tem- 


peratures), recognizes that quantum spins of spin quantum number S can orient in 2S+ 1 


directions in real space. The spin vectors have length 4/5 (S +1) in units off. We now 
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replace the quantum spins S; > VS(S +1)s, by classical vectors of length JS(S ely, 


where §; iS a unit vector at site 7, free to point in any direction. 


The classical Heisenberg Hamiltonian is now 


N N 
¥ =-JS(S+1)>°s, -8,,, —mBJS(S +1)> s? . (3.2) 
el i=] 


To simplify the notation, establish effective dimensionless coupling constants and 


create simplified Boltzmann factors, we will define K = B/S(S +1) = JS(S +1)/(k, 7), 
where J is the absolute temperature and k, =1.38x10J-K~'. Similarly, we define 


L = BmBJS(S +1) = mB./S(S +)/k,T . The Boltzmann factor, from statistical 


mechanics is then expressed as, 


exp(- Bx) = exo KD. Ss ix =| (3.3) 


i=] | 
The partition function as the sum of the Boltzmann factors becomes an integral 


since the classical spins are continuous, 


Z = || — [T]4s,exp(- x)= one {Tas oa KD Sa +t¥s1), (3.4) 


S, Sy Fl S; Sy Fl 
where ds, = sin0,d6,d@, , is the element of solid angle about “spin” s;. In parallel with 


our treatment of the Ising model in Chapter I, since the spins are classical variables, and 


hence commute with each other, we can express the exponential of the sum in (3.4) as a 


product of exponentials, 
; 7 
Z= | a Nalas;2 Gass ih Gass) (curse \(sy58,), (3.5) 
S, Sy #! 
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where 7 (s, s’) is the transfer operator. In what follows, we examine some of the 


mathematical properties of the transfer operator. We then show that knowledge of the 


eigenvalues of 7 enable us to find the partition function. 


One way to wnte 7 (s, s‘) is aS a Symmetric function of its arguments, 
, 4 ] 76 92 
T(s,s’)= exp(Ks-s +s +S }. (3.6) 


This form is not necessary, however; one could also contrive the following non- 


symmetric transfer operators: 
T. (s,s’)= exp(Ks -s'+ Ls* | or T(s,s’)= exp(Ks -s'+ Ls'*): (sg) 
these non-symmetric operators are in fact transpose pairs since 7, (s’,s)= T_(s,s’). 


In Chapter IIL, Section E we will return to these non-symmetric forms of 7. Here we will 
explore the mathematical consequences of employing the symmetric version of 7, (3.6). 


If we consider the integral equation 
Jas'T6s,s’Wv,(s') = Awl), (3.8) 
this defines the eigenfunctions y,(s) and eigenvalues A, of 7, where alo cee re 1S a 
discrete index. If 7(s,s’) is symmetric, (e.g. equation (3.6)), then Hilbert-Schmidt 
theory [Ref.29] guarantees that the eigenvectors are a complete orthonormal set and that 


the eigenvalues are real. The completeness and orthonormality of the eigenvectors on the 


unit sphere means: 


[as y*(s)w,,(s)=6,,, (orthonormal) and (3.9a) 


Svi(s) y.(0')=6(s~ 2’) (complete) (3.9b) 


n=] 


It is shown in the theory of integral equations [Ref.29] that because (3.6) is real and 
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symmetric, we may expand the transfer operator in terms of its eigenvalues and 


eigenvectors, 
T(s,s')= DA, (syn (s'). (3.10) 
n=] 


With this expansion, one can show that the partition function for the N-spin system is 


given by, 
ae 3.11) 
n=) 
To see how (3.11) arises, let us work out the details for the case of the N=2 


system. Starting from (3.5) we have 


Z, =| \idsas, 15,5, \0(s,5, )). (3.12) 


5, 52 


Substituting the expansion (3.10) in (3.12), we have 


Z, 


J fads, © Av 62 (62) AW a(S Wals) 


5} 52 


y ham [ 4s, Y (s, Ws (s, fds. (s, WV, (s, )= 3 44m nm 


n or=) n,m=) 


Z,=> A,, (3.13) 


where we have used the orthonormality properties given in (3.9a). It is thus clear how to 
extend the treatment to general values of N and arrive at (3.11). Just as we obtained in 
the analysis of the Ising model, obtaining the partition function of the classical 


Heisenberg system 1s tantamount to finding the eigenvalues of the transfer operator. 


va 


33 


C} ZERO-FIELD EIGENVALUES 


At this point, it will be useful to show how one may obtain analytically the 
eigenfunctions and eigenvalues of 7 for zero applied magnetic field. The key to the 


subsequent development is the use of the following expansion 


exp(Ks, -s,)= 42, Xf K Fon (5,)¥7, (52), (3.14) 
(a 
where f(K)= ag 1uy(®), (3.15) 


is a modified spherical Bessel function and Y, ,, is the standard spherical harmonic 
function. For future reference in characterizing matrix symmetry, we note the Bessel 
function parity property f,-K)=(-1) 7, (K ). In addition, another essential result, 
commonly called the Condon Shortley phase convention, is the fact that the spherical 
harmonics obey Y,,, = (— 1)” Y, _,, - In what follows, we will work with spherical polar 
coordinates (0, p). The angle between the two “spins” (unit vectors) 1s given by 

S,-S, =cos@ = cos@, cos6, +sin8@, sin 8, cos(¢, — ¢, ), (3.16) 
which is a standard result from vector analysis. We will use the shorthand notation 


iG (s) to denote a (0, ¢). Equation (3.14) follows from combining two results from 


mathematics. The first is the expansion, 


exp(z cos 6) = x (27 +1)f,(z)P (cos), (3.17) 


which is a special case of the Gegenbauer addition theorem for Bessel functions. 
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The second is the addition theorem for spherical harmonics, 
P(cos®) = oY 6,6) ¥, @2,92). (3.18) 
m=—!| 


Combining (3.17) and (3.18) leads to (3.14). 


Comparing (3.14) with (3.10), we can identify the eigenvalues of T as 


Atm =4af; (K). Note that in this zero-field case each eigenvalue is (2/+1)-fold 


degenerate. When we turn on a B field, the degeneracy is lifted. We also identify the 


eigenfunction of T in the zero-field limit as y, ,, (s) =Y!,,, (s). Note that, as opposed to 


(3.8) in which a single generic index is used to label the eigenfunctions, in this case we 
must employ two indices to label the eigenfunctions and eigenvalues. We can verify that 


the Y, ,,’s are the correct eigenfunctions with the /, (K)’s as the eigenvalues as follows. 
Let’s assume this assertion to be true and substitute into the eigenvalue integral equation, 
J exp(Ks, «82 )¥ pp (S_ 5. = Ap Vp (S, ). | (3.19) 
Now employing the expansion (3.14), and the (3.9a) orthonormality property we have 
4n| ds, DSK Wm (5 Fin 2 Berm 82)= 4 LK Fin ) 452% 0 (S2 Vern (52) 
= 405, Ain Bran = 4h Kral) _— 
Thus, as Fenonate by Joyce [Ref. 7], in zero applied field, the partition function is 


Z,(K)=(42)" > @+1) AN (K). (3.21) 


We can check that (3.21) properly reduces to the correct high-temperature limit 


Z(0)= (4) upon using the property of the modified spherical Bessel function 


S (0) =010- 
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D. BLUME, HELLER, AND LURIE METHODOLOGY 


Published in 1975 [Ref. 8], this numerical method achieved the first theoretical 
results of the thermodynamic properties for the classical Heisenberg magnetic chain in an 
applied magnetic field. Using the transfer matrix method and numerical Gaussian 
integration, this achievment extended the zero field analytic solution of Joyce to the non- 
zero field case. In a magnetic field, formulating a numerical evaluation of the 


eigenfunction, Blume et al [Ref. 8] first noted that the eigenfunction can, by symmetry, 
be written in the form ¥, _(s)=y,,, (cos Y | , where @ and ¢ are the polar 


and azimuthal angles respectively of spins. This separation of variables recognizes that 
even in the presence of an external field, the azimuthal parameter m remains a good 
quantum number. Letting x =cos@, utilizing (3.16), and using the symmetric kernel 


(3.6), equation (3.8) (with indices /m in lieu of 7) becomes, 
1 Tas. 3 1 
[ dx’ [ 2g’ exp (Kes ¢ l( _ x? \ - x’? )p cos (6 — ¢’)+ = (x + «| 
-1 0 : 
(3.22) 


imo’ 1 
uae «= een (=)ome'*? 


The integral over ¢’ can be carried out analytically with the result, 
1 IE 1 

21 | ax’ exp{ Ki’ + 5 + NK li — x? \i — x” MV Jone) =A Wim(x). (3.23) 
=I 


Le 
Here J, (x)= = exp(xcos¢ —im¢)d¢d =I,,(-x)=L,, (x) is the Bessel function of 
Md 0 


imaginary argument. This one-dimensional integral equation can now be solved 


numerically by converting it to a matrix equation. The integral over x’ is performed by 
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N 
.: f 
N; -point Gaussian integration, using the approximation [, f (x}& = >» We (x.), where 
5 =; 


the weights w; and the points x; are tabulated.[Ref. 30] The integral equation then 


becomes, 


N 
[_&'G,, (x, x" W in (x') = Ain im (x) = >'w.G, (x, x, Wp (x, ) (3.24) 
fe. 


ae 
where G,, (x, x’)= ne 1 = (x ( — x" fl —x”? \}? . If we look for solutions of Eq. 


(3.23) only at the points x = x, of the numerical integration, the integral equation 


becomes a matrix eigenvalue equation: > w Ay (x 25 VW (x = Ain Win (X;). To make 
j 


this more symmetric, we multiply both sides of the equation by JY; , obtaining 
ay 


DON” = Ani” » (3.25) 


j=l 


where Ae = Jw, G,, (x,,x, fw, and gi) — JY; Win (x. ). Equation (3.25) is an 
N,xN, matrix eigenvalue equation with N; determined by the number of points used in 
the numerical integration. As stated in the introduction, (and assuredly worth repeating), 
a value of N;=/6 suffices to give convergence to seven significant figures for all values 
of K and L. The largest eigenvalue of equation (3.25) occurs for m=0. This is in 
consonance with the fact that the eigenfunction belonging to the largest eigenvalues has 
no nodes. In a sufficiently large spin system, only the largest eigenvalue survives, that is 
Z-—A™ for N-»>« . The free energy appropriate to an infinite-site system can be 
obtained from this eigenvalue and other thermodynamic quantities can be found by 


numerical differentiation with respect to the appropriate variables. 


sf 


E. REPRESENTATION IN SPHERICAL HARMONICS 

Mark Auslender, an Israeli physicist, has suggested, in a private communication, 
[Ref. 9], an alternate numerical strategy for solving the eigenvalue problem associated 
with the transfer-operator, (3.8). Auslender’s suggestion consists of representing the 
transfer operator in terms of a spherical harmonic basis set. Since the spherical 
harmonics are the eigenvectors of the transfer operator for zero magnetic field, it is 
possible that Auslender’s suggestion could prove numerically more efficient or flexible 
than the Gaussian integration approach discussed above for the case of non-zero magnetic 
field. We stress that Auslender has merely suggested this approach; we are here working 
out the details. 


Thus, we first expand the eigenfunctions of 7 in terms of spherical harmonics, 


VY, (s)= DOr Yom (s). (3.26) 
Substituting (3.26) into (3.8), we have 


C1, fas 76,5 Fru) =A, DOT a) 3.27) 


Then, multiplying (3.27) by Y,,, (s) and integrating over s, we obtain 


DCm [as[ as’ 7”, (s)T(s, se (sa) a "ean (3.28) 


We now define the matrix elements, 
T(, mM, is m’) is [ [asas'y,”, (s)r(s, OV (s’), (3 .29) 
in terms of which (3.28) becomes 


Dea! myen A.C". (3.30) 
I'm’ 
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For each m=1,2....., (3.30) has the form of a matrix eigenvalue problem. In what follows, 
we will simplify our notation and suppress the index n. As previously, the goal will be to 
obtain the eigenvalues of 7, now in a matrix representation given by (3.29). 

The form of (3.30), however, is slightly unusual because the components of the 
eigenvectors are labeled by a double index set. We can cast (3.30) into standard form by 


lining up all of the eigenvector components into a column vector, 
C= “0 | (3.31) 
C 


The column vector C is infinite dimensional and the eigenvalue condition becomes a 
matrix eigenvalue equation, T-C = AC, where T is an infinite-dimensional matrix, the 
first few elements of which would be, 


T(0,0;0,0) 7(0,0;1,1) 7(0,0;1,0) 7(0,0;1-1) 7(0,0;2,2) --- 
T(1,1;0,0) 7(11;1,1) 7(1;,0) 7(,1;,-1) 7(11,2,2) -- 


T(1,0;0,0) 7(1,0;1,1) 7(,0;1,0) 7(1,0;1-1) 7(1,0;2,2) --- 
‘les (3.32) 
RC 1,00) 0714.1: 1 ieee = 1) (122) oe 


eis (APRN! JE AIeUO)| TP Pee) 1090-22) 


We have thus formally reduced the integral eigenvalue equation (3.8) to a matrix 
eigenvalue equation,where, however, the matrix is infinite dimensional. Clearly, some 
practical means of truncating the matrix must be developed. This remains to be done. 


We note that, in principle, the Gaussian integration scheme also produces an infinite 
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matrix to diagonalize. There it was found empirically, that a 16 x 16 matrix produced 


satisfactory results. 


It will be instructive to first evaluate the matrix elements of T for the case of zero 
magnetic field by substituting (3.14) in (3.29) 
T(,m,I',m')= | ds | ds'Y,”,(5)T(s, 5’ Yn (S’) 
= | ds| ds'Y, (s)exp(Ks-s’)¥,.,(s’) 
= ‘/ Fi, (K)| as | ds"Y,7, (S)Y, mm (SWiiny (S’ Worm (5’) (3.33) 
=47, ef, (K)OmOnmO; One 


hm, 


= An f, (K)6, » Overt: 


Not surprisingly, yet reassuringly, this shows the matrix is diagonal with the zero field 
eigenvalues along the diagonal. 

To incorporate the magnetic field into the matrix construction, either a symmetric 
or a non-symmetric approach is possible. We will explore both avenues, but before 
proceeding, we list the following results that will prove useful in our subsequent analysis. 


The first derives from the addition theorem for Legendre polynomials: 
exp(Ls‘) = exp(Zcos@) = V4z > J2j +1 ve UO (s). [Ref 30 p.445] (3.34) 
j=9 


We will also require the general integral over three spherical harmonics, 


(22 EY ORY ORT CL a es 
dsY uy pa le 2 emcall(as 
J a (5) fava (Ym (s) Ar t 0 af m, al ) 


aa ] 
where j : is the Wigner 37-symbol. [Ref. 31 p63] The Wigner 37 symbol 
m Mm, my) — 


is a symmetric form of the Clebsch-Gordan coefficient that arises frequently in contexts 
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involving coupling of angular momenta in quantum mechanics and other applications of 
the rotation group. The 3j symbols are non-zero only when m, +m, +m, =0 and the top 
row Satisfies a triangle condition, |/, -/,|<J, </, +1,. When the bottom row is 
identically zero, there is an additional rule that the top row must sum to an even integer. 
1. The Symmetric Kernel 
Since (3.6), the original transfer operator 7 (sus) is symmetric, we expect the 


matrix T to be Hermitian. This will ensure that the eigenvalues A are real. Following 


(3.29), we need to evaluate the matrix elements, 


T(,m:Um!) = [asf as'Y (9I0(5,8 Wen) 


3.36 
= | ds | asia (s)exp(Ks -S'+ LL (s? +s" Wa (s’). ade 
We now substitute in the expansions (3.14) and (3.34), 
exp(Ks-s’)= any FAKY,,,(s)¥,”,(s’) and | (3.14) 
I,m 
exp(L Ls?) = Vaz > f,(L/2W2/+ IY, .(s). (3.34) 


1=0 


The expression for the matrix element then becomes, 
(1, m;1',m’)= | ds{ ds’Y,”,(s)exp(Ks-s'+4L(s? +5’? 4 (s’) 


= fasfasy (faa fl (Wm 0 
¢ RS ALIANT, )] ARTI, CAME tl) 0. 


1,=0 


4] 


Collecting terms, we have, 


T(L,m,I',m')= (x) 2d Ff, (K)f, (L/2)Y21, +14, (L/2)/ i, +1) 
hem |, ds E87) 


x ([as¥ 0 (Mam Ono) fa5'F.n 6,06 rls’). 
Employing (3.25), the integral over 3 spherical harmonics and using the fact that 


Y, = (-1)”¥,_,,, we obtain the ponderous expression, 


T (1, m;I',m’)= 3p ap apd Ff, (K)A, (L/2)Y 2h, +f, (L/2)Y(22, +1)” C1)” 


Am i, J, 


(Sane  +:1K22, My" CL TIE GE I 
4n 00 0;\-mm, 0 
fee eer +a)" hi,Vy¥h &P? 
An 00 Ol\-m 0 m') 
We can simplify this beast somewhat, 


T(,m;I',m')= (42) >> f, (KDA, (L/2)A, (L/2 X22, + 1X20, + 1X20, + 1)Y (27 + 1X2!’ +1) 


I mh L od, 
’ l le ies l bs Ne ; : Me i l L (- ne 
00 O\-m m 0)\0 0 Of-m Om 


Utilizing the properties of the 37 symbol discussed above, we have for the symbols 


(3.38) 


appearing in (3.38), 


We can thus instantly sum over mj, (using the fact that 2° = 6,,,"), and the 


m, m7 Ona 


42 


matrix elements become, 


T (i, m;I',m’) = 476,__,|(21+ VCS) > y FAK) A, (L/2)f, (L/2X20, + 1X2/, + 1X22, +1) 


ik 1,=0 1,=0 
‘alale Commas es alee ‘ome! 5 r aaa e Er 
x 
0 0 Ofh-m m 0 0 O Oh-m O Mm 


Equation (3.39) seemingly involves a triple infinite sum. The “triangle” properties of the 


(3.39) 


37 symbols probably restrict some or all of these summations to be finite sums. We will 
not further analyze (3.39) because, as it turns out, a considerably simpler expression for 
the matrix elements can be obtained by working with the asymmetric form of the transfer 


operator (3.7) mentioned in chapter IIL section A. 


2 The Non-Symmetric Kernel 


We now obtain the matrix elements associated with the asymmetric form of the 
transfer operator, 7(s,s’)= exp (Ks -S'+ Ls”). First, however, we note that one might be 
concerned that a non-symmetric kernel would not have real eigenvalues. (Recall, as 
discussed above, that Hilbert-Schmidt theory guarantees that a real symmetric kernel has 
real eigenvalues). In this particular case, however, we can show that the non-symmetric 
kernel is related to the symmetric form of the kernel by a similarity transformation, and 
hence has the same eigenvalues as the symmetric kernel. Consider that the eigenvalue 


equation, (3.8), is equivalent to the following 


[s(.s')e7 sels’, ()as’ = 2,.2(sv, (5) (3.8) 
or le (s,s’)¥, (s’\ds’ =A VY, (s) 


where T(s,s’)= g(s)T(s,s‘)g"(s') and ¥,(s)= (s)y,(s). 
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So long as a function g(s) can be found such that T is symmetric, Zand 7 will have the 
same (real) eigenvalues. In our case, it is easy to find the transformation function, 
g(s)= exp Ls). Thus, the (3.7) non-symmetric kernels 7, (s,s’) have the same 
eigenvalues ” the (3.6) symmetric 7(s,s’). 
Utilizing now T, (s,s’) from (3.7), we have from (3.29), 

T(l,m;1',m')= | [dsds'Y,”,(s)exp(Ks-s’ + Ls? Vy,.,.(s 

= (2 OSA, sain (s\n, (Y,.0(8) 

«(fas'ys,, (5° Wrw(s’) (3.40) 
= (tn) 5 (R)S I, (ON2F A Yn OY, 0) 


= sf NOTA DIS 4s iY 4 4 le er 


00 0 mim’ 0 
where, in arriving at the last line, we have used the phase convention that 
Y= (—1)” Y,_,- We now utilize the “selection rule” properties of the 3j-symbols: The 
3j-symbols are non-zero only when the sum of the bottom row 1s zero, and when the 
upper row Satisfies the triangle inequality, H —] i < j <1+1'. These two facts: (1) restrict 


the sum over j,and, (2), makes the matrix element diagonal in the variable m: 


mn Pegy iP j 
T(L,ml',m')=42\-1) 6 K )f (27 +127’ + LYN2j7 +1 
(sm m= AeA 5, ARNON FC) 5 fo 8] 
(3.41) 
We note the additional rule for 37-symbols that when the bottom row 1s all zero, the 


symbol is non-zero only when the upper row sums to an even integer. Thus, the sum 


over j in (3.41) is further restricted to values such that 7+/'+ j =2p, where p is an 


integer. We note that the fact that T(J, m;/',m’) is diagonal in m, makes sense; the 
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azimuthal quantum number remains a good quantum number in the presence of the field. 


Finally, we note that utilizing the other asymmetric form for the kernel 
(rT (s,s’) = exp(Ks -s' + Ls’? ) merely produces the transpose of T(J, m;1',m') 
It is easy to verify that the correct zero-field limit results from (3.41). Using the 


fact that f, (0)=6 jo and the 3/-symbol, 





( 1] 4 _@I” 
-mmo0) ial 
we obtain the zero-field limit, T(J, m;/ 'm') = 47f, (K So as derived previously (3.33). 
One can also show that (3.39) properly reduces to (3.33) in the zero-field limit. 
Using the symmetry properties of the 3j-symbols, it is simple to show that the 

transpose of 7 (i, m;',m') is given by, T(i',m',1,m)=[f,(K)/ f, (K) (I,m:1',m'). 

Since we started from an asymmetnic version of the transfer operator, it is not surprising 
that T is not symmetric. We note that the transpose relation implies that the eaten 


elements have the symmetry property, 7,7 (/ "ms, m)= FT (1, m;l m’') 
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IV. DISCUSSION 


The purpose of this thesis has been to investigate a method for obtaining the 
thermodynamic partition function for classical Heisenberg spins that interact with 
isotropic nearest-neighbor exchange interactions and which are coupled to an external 
magnetic field. Equation (3.41) is the main result of this thesis. It provides an expression 
for the matrix elements of the transfer operator for classical Heisenberg spins that result 
upon utilizing a basis of spherical harmonic functions. The motivation for pursuing this 
new matrix representation is ultimately to assess its numerical efficiency, as compared 
with Blume’s Gaussian quadrature method, in determining the eigenvalues of the transfer 
operator. Within the transfer matrix formalism for calculating the equilibrium properties 
of interacting spins on a lattice, the “matrix” in this case being an operator, the partition 
function is obtained from the eigenvalues of the transfer operator. We emphasize that the 
Gaussian integration method 1s the only other numerical method available in the physics 
literature of which we are aware for obtaining the eigenvalues of the transfer operator for 
classical Heisenberg spins in an external magnetic field. Moreover, we note that the 
transfer matrix method 1s the only general method for treating the statistical mechanics of 
interacting spins; there are only a handful of exceptional cases where the partition 
function can be obtained directly, without recourse to the transfer-matrix method. 

As discussed in the introduction, great progress is being made in the ability to 
fabricate molecular clusters containing a small number of magnetic ions (e.g., as few as 
four). Recent advances in the fabrication of molecular magnets portend an 
unprecedented ability to control the placement of magnetic moments in molecular 


structures and hence to design and produce nanometer-scale magnetic systems. As 
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molecular magnetic systems continue to be explored for their possible applications, 
robust and reliable numerical methods will be required to model their thermodynamic 
properties. In the past, physicists have explored spin models for their ability to 
characterize phase transitions. Phase transitions and critical phenomena, however, require 
that the “thermodynamic limit” be taken at the end of the calculation, which, for spins on 
a lattice, means that the number of spins in the system becomes infinite, (i.e., VW > o ). In 
this limit, only the largest eigenvalue of the transfer matrix becomes relevant. We note 
that specialized numerical methods exist for seeking either the smallest or the largest 
eigenvalue of a given matrix. For the development of nanomagnetism, however, we are 
concerned with the opposite limit to that attendant to the study of phase transitions, (1.e., 
here N-> finite few). To obtain the partition function for systems with just a few magnetic 
atoms, we will require an indefinite number of the transfer matrix eigenvalues and hence 
it is worthwhile to explore new methods for calculating these quantities. We note that the 
number of eigenvalues of the transfer operator 1s independent of the size of the system. 
In some sense, we have entered an era of “applied statistical physics,” and appropriate 
tools are required. 

Without a detailed numerical investigation, it is difficult to assess the utility of 
(3.41) vis-a-vis Gaussian quadrature. We can offer the following observations. First, it is 
exact. Equation (3.41) provides the exact matrix elements of the transfer operator in the 
basis of spherical harmonics, and this fact alone may offer insights. Stated differently, 
the Gaussian integration method is purely numerical, whereas (3.41) is based on an exact 
theoretical expression, which in and of itself may prove useful. Second, one hopes that 


(3.41) will prove advantageous at least for the case of relatively small magnetic fields. 


48 


Because (3.41) is diagonal for zero magnetic field, (since the spherical harmonics are the 
zero-field eigenfunctions of the transfer operator), one would expect for non-zero 
magnetic field that the off-diagonal terms would remain relatively small, and, moreover, 
to become progressively smaller as one proceeds away from the diagonal. This follows 
from the properties of the modified spherical Bessel functions in (3.41). These functions 
have the property that they become monotonically smaller as a function of the order for 
fixed values of the argument. In particular, when the value of the order exceeds the value 
of the argument, the value of the function vanishes (approximately) exponentially as an 


increasing function of the order. From (3.41), we see that the lower limit of the 


, 1.e., the order of the first (and largest) term in the 





summation is given by |/-/' 
summation 1s directly given by the distance to diagonal. Hopefully, such considerations 
will prove useful in developing “rules of thumb” for deciding how to truncate the matrix 
for the purpose of numerically obtaining the eigenvalues. In a similar way, we note that 
the number of terms to include a/ong the diagonal is governed by the overall, field- 
independent, modified spherical Bessel function in (3.41) which is a function of K, the 
dimensionless nearest-neighbor coupling constant. It thus seems likely that with suitable 
numerical experimentation, one can develop practical schemes for truncating the matrix 
for given values of K and LZ. Finally, a decided advantage of (3.41) is that it provides a 
systematic way for increasing the accuracy of the eigenvalues, if such is desired. It was 
noted in Chapter II that a 16x16 matnx was sufficient to guarantee seven digit accuracy 
with Gaussian integration. It is a feature of Gaussian integration, however, that one 
cannot systematically obtain more accuracy by increasing the number of integration 


points. For numerical reasons, the accuracy of Gaussian integration “saturates” for a 
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relatively small number of integration points. Thus, while seven-digit accuracy is 
commendable, if for some reason one wanted higher accuracy, it could probably not be 
obtained using Gaussian integration. 

As surmised at the outset, only performing the actual calculations will convey 
“the rest of the story”. In general, computational efficiency yardsticks, like the means of 
matrix truncation and the ends of convergence results, must of course await actual 
programming and calculation. Application accuracy, as well as flexibility are key factors 
in adjudging the utility of any tool. So it is with models; in this case the efficacy, 
versatility, and robustness of the spherical harmonic representation is yet to be 
determined. Although pragmatic results will remain the preeminent objective, when an 
intrinsically exact and high fidelity model formulation can engender an illustrative 
understanding of the phenomenon examined, this is a welcome bonus. Finally, as for the 
macromolecular magnetic frontier, the nanometer investigation and fabrication 
technologies undoubtedly will encourage a convergence of scientific disciplines. The 
chemist, physicist, and “nanotechnologist” will merge here to both gain an understanding 


and inevitably utilize these ultra-small complexes. 
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